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We study the quantum entanglement caused by unitary operators that 
have classical limits that can range from the near integrable to the completely 
chaotic. Entanglement in the eigenstates and time-evolving arbitrary states 
is studied through the von Neumann entropy of the reduced density matrices. 
We demonstrate that classical chaos can lead to substantially enhanced en- 
tanglement. Conversely, entanglement provides a novel and useful character- 
ization of quantum states in higher dimensional chaotic or complex systems. 
Information about eigenfunction localization is stored in a graded manner in 
the Schmidt vectors, and the principal Schmidt vectors can be scarred by the 
projections of classical periodic orbits onto subspaces. The eigenvalues of the 
reduced density matrices are sensitive to the degree of wavefunction localiza- 
tion, and are roughly exponentially arranged. We also point out the analogy 
with decoherence, as reduced density matrices corresponding to subsystems 
of fully chaotic systems are diagonally dominant. 
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I. INTRODUCTION 

Entanglement has been studied since the early days of quantum mechanics and provides 
the essential ingredient of EPR phenomena Jl]. It results when the state of a system, which 
is composed of at least two subsystems, cannot be written as a product of states that reside 
entirely in the subsystems. This leads to the well known unique quantum correlations that 
exist even in spatially well separated pairs of particles. More recently entanglement has been 
discussed extensively in the context of quantum information theory. It has been recognized 
as a quantum resource for quantum superdense coding and quantum teleportation, while 
helping make quantum computations qualitatively superior to classical ones 0. 

Usually the entangled subsytems are distinct identical particles and the entanglement is 
created by symmetrization, such as in the spin singlet state of a fermionic pair. However, 
entanglement may of course be created during time evolution due to conventional interactions 
of a potential nature. In this case, the "subsystems" may also be the many degrees of freedom 
of the same particle. This is more general than the special case wherein the interaction 
preserves the permutation symmetry between the say d (formal) degrees of freedom and 
may therefore be interpreted as representing d identical particles in one dimension. 

Presented with the inexpensive resource of unentangled states, unitary operations may 
be devised to create potentially useful entangled states. Thus an understanding of the range 
of entanglement produced by the wide range of unitary operations is of interest. If these 
unitary operators are generated from Hamiltonians, then a natural question that arises is 
the connection between the dynamics generated by the Hamiltonian and the entanglement 



produced. In particular it has been recognized for some time now that the two extreme cases 
of classical dynamics, namely the completely integrable and the completely chaotic ||, leave 
remarkably different traces on quantization In this paper we enquire into the possible 
effects of quantum chaos on quantum entanglement. A related issue is that of decoherence 
H and our treatment is equivalent to treating the environment as a subsystem of a quantum 
chaotic system. 

Certain finite unitary matrices have been studied extensively in the quantum chaos litera- 
ture and represent a rich source of operators with a wide range of well understood dynamical 
behaviors in the classical limit .They may for instance arise naturally while quantizing sym- 
plectic maps on the torus 0-§|. However, most of the studies have been limited to the 
two-torus, with essentially one-degree of freedom, while entanglement studies necessarily 
deals with the entanglement between two subsystems. Thus, while two-degree of freedom 
Hamiltonian flows lead to one-degree of freedom maps, we need to consider at least a four- 
dimensional symplectic map to study entanglement directly. 

We use a coupled quantum standard map as the model system. The standard map is 
well understood classically, there is a smooth transition from the regular to the chaotic, and 
has been the object of several quantum investigations as well (reviewed in ||), including 
an experimental realization M. The coupled map has been studied classically, principally, 



to understand Arnol'd diffusion |]i0| . Thus our results below also represent new results 
in the largely sparse literature on higher dimensional (d > 1 for maps, d > 2 for flows) 
quantum chaotic systems. In particular we explore the usefulness of entanglement as a 
measure of quantum chaos for such systems. Our results, largely exploratory and numerical 
in nature, quantify what is intuitively (and certainly linguistically) evident: the connection 
between entanglement and chaos. The larger the chaos, more the entanglement that is 
induced by the dynamics, as opposed to entanglement due to symmetries. However we also 
find interestingly that the reduced density matrices store information about localization of 
states in a systematic manner. 

Consider a bipartite quantum system whose state space is H = TCi <g> H2 where 
dim7ij = N, (i = 1,2). The unitary operator acts on vectors in the space 7i where 
dim7i = N 2 . We consider the case where the two subsystems have the same dimensional- 
ity. Let \ip) eH, and the reduced density matrices be p\ = tr 2 (| , 0}(V'l) an d Pi = tri(|-0) (ip\), 
where the first matrix is obtained by tracing out the second degree of freedom and the second 
by tracing out the first. The von Neumann entropy S, referred to in the rest of the paper 
as simply entropy, of the reduced density matrices measures the entanglement of the pure 
state in an essentially unique manner. 

S = -tri(pi log(pi)) = -tr 2 (p 2 log(p 2 ))- (1) 
The close analogy between this entropy and thermodynamic entropy has been noted and 



discussed earlier [11 



The eigenvalues of the reduced density matrices form the (square of the) coefficients of 
the Schmidt decomposition. This decomposition expresses the state in the full Hilbert space 
as a linear combination of N product states instead of the evident N 2 . If we were to write 
the eigenvector for instance in the position basis and in the Schmidt decomposed forms we 
would have 
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N 

m = E w>i>M = £V^lfc>Vi> 5 - (2) 

m,ri2 j=l 

A simple proof of this is found in the Appendix of |L2] . Thus there is a basis in each Hilbert 
space such that each basis vector in Hi is uniquely correlated with one in 7i 2 as far as the 
state \ip) is concerned. The states \4>j) and \4>'j) s are the eigenvectors of the reduced density 
matrices, and form the basis of the Schmidt decomposition. The Schmidt decomposition is 
not merely a mathematical convenience but is understood to provide a deeper understanding 
of correlations between subsystems. 

In the literature on quantum chaos or localization H,[yJ , the Shannon entropy has often 
been calculated, and for the above state it is given by: 

S s han(\lp)) = ~ £ l C nm 2 | 2 lo g(l c nin 2 | 2 )- (3) 

ni,n 2 

This quantity is of course basis dependent and vanishes if the basis is chosen to have a direc- 
tion along The von Neumann entropy or entanglement is also basis dependent, but is to 
an important extent immune to arbitrariness by being invariant under local unitary trans- 
formations, even such not being the case for the Shannon entropy. Local transformations 
are those that act only on individual subspaces. 



II. FOUR DIMENSIONAL STANDARD MAP 



We now define the four-dimensional standard map JTIJ. It is composed of two pendulums 
that are periodically kicked and are also coupled to each other. The symplectic transforma- 
tion of the phase space variables (q%, q2,Pi,P2), connecting states just before two consecutive 
kicks, separated by unit time, is the classical map: 

q'l = 9i + P[ (4a) 

K b 
Pi =Pi + 7T sin(27rgi) + — sin(27r(gi + q 2 )) (4b) 

<?2 = <?2 + P' 2 (4c) 

K b 
P ' 2 = P2 + 2tt sin ( 27T,?2 ) + 2^ sin ( 27r (<?i + <?2))- (4d) 



The phase space is restricted to the unit four-torus T 4 , and therefore mod 1 operations 
are understood in all of Eqs. (Q). If b = 0, the system falls into two uncoupled standard 
maps. In this limit much is known of the dynamics 0; briefly, if K = (referring now to 
either K\ or K 2 ), the dynamics is integrable, while at K w 1 the last KAM rotational tori 
breaks, heralding large scale diffusion in phase space, for K < 5, the phase space is that of 
a typical Hamiltonian system, a mixed phase space with both regular and chaotic regions. 
When K » 5, the dynamics is practically completely chaotic with possible appearances 
of very tiny stable islands through tangent bifurcations. For 6^0, little is known, due 
to the dimensionality of the phase space, although this map has been used in studies of 
Arnol'd diffusion. We suggest, and substantiate below, that in cases such as these where 
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finite unitary matrices may be constructed as quantization, the quantum maps can be used 
to actually find transitions to classical chaos. As we noted for the general case earlier, in 
the case K\ = K 2 , the system possesses permutation symmetry between the two degrees of 
freedom and may be interpreted as two interacting particles in a one-dimensional standard 
map external potential. 

The quantization of the symplectic transformation in Eq. (f|) is a finite unitary matrix, 
whose dimensionality is A 2 , and A = 1/h, where h is a scaled Planck constant. The 
classical limit is the large A limit. The quantization is straightforward as there exists a 
kicked Hamiltonian generating the classical map. The quantum standard map on the two- 
torus in the position representation is 

U{n',n;K u a,P) = exp [~ iN -^ cos (j^( n + «))) 

x ]T exp + P) 2 ) exp (^(™ + P){n- n')) . (5) 

The position kets are labeled by n = 0,N — 1 and the position eigenvalues are (n + a)/N 
while the momentum eigenvalues are (m + /3)/N, m = 0, . . . , N — 1. Here a and j3 are real 
numbers in [0, 1) which represent quantum boundary conditions and are convenient devices 
for breaking phase space reflection symmetry (the phase a) and time reversal symmetry (the 
phase f3). The four-dimensional quantum map is but a simple extension: 

(njTialWlniTk) = Uirf^n^ K u a, (3) U(n' 2 ,n 2 ; K 2 , a, 0) 

x exp f -iN^- cos(^(ni + n 2 + 2a)) ] . (6) 
y Z7r JM J 

IA is a unitary matrix in T~C, and will induce mixing between the two subsystems. We have 
assumed the quantum phases in both the subsystems to be identical. Throughout this paper 
we use a = 0.35, (3 = as the quantum phases. 



III. RESULTS 

A. Stationary State Properties 

If b 7^ an unentangled initial state, such as \rii,n 2 ) = |ni) <S> \n 2 ), would eventually 
get entangled by the repeated action of the unitary operator U. The entangling proper- 
ties depend on the entanglement already inherent in the stationary states or eigenstates 
l^i), i — 1, . . . , A 2 , of U. Thus we first calculate the average entropy of the eigenstates 
when Ki = 0.1, K 2 = 0.15 as a function of the coupling constant b. At these values of Aj 
the uncoupled standard maps are almost wholly regular. We calculate 

_^ N 2 

S = tT i(pu l °g(pii))i Pu = tr 2 (|^i)<V>i|) ( 7 ) 

iv i=l 

where pu is the reduced density matrix after tracing out the second degree of freedom. The 
quantity S is a gross quantity averaged over the entire spectrum, and gives an idea of the 
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average entanglement we can expect on using the operator U. In Fig. (|]) we see the entropy 
increasing from zero at 6 = and attaining a nearly constant value beyond 6 ~ 3. 

The increase in the entropy proceeds along with a gradual increase of chaos in the system, 
flattening out after considerably uniform chaos has been achieved, a fact that is confirmed 
by iterating the classical map Eq. ([|). In fact this suggests a compact way of exploring 
the classical transition to chaos which is otherwise mired in problems of visualizing four- 
dimensional sections. Thus entanglement is clearly a function of the nature of the underlying 
classical dynamics. That the entropy increases as an approximate power law before flattening 
out is shown in the inset. Roughly we get S ~ b A . We note here that the Shannon entropy 
would behave differently, as even when 6 = there is a non-zero Shannon entropy in general. 
If there is large scale chaos in the subsystems (as in our model if K\ and K2 are greater 
than five) it will be reflected as a large Shannon entropy; however the entanglement would 
be zero. For large coupling between the subsystems, both the entropies appear to be well 
correlated. 

Subject to the constraint that tri(pij) = 1, the maximum entropy is log(iV), and corre- 
sponds to the "microcanonical ensemble" with all the eigenvalues (Xj, j = 1, . . . , N) of pu 
being equal to 1/N. The entanglement entropy induced by the dynamics of quantum chaos 
falls short of this and in fact the eigenvalues if arranged in decreasing order are exponential 
and reminiscent of the "canonical" ensemble; see discussion below. From the data of Fig. ([!]) 
it appears that at saturation S ~ log(0.59iV). Thus roughly 0.59iV pairs of correlated states 
from the two subspaces make up a typical state. The meaning of this is that of the optimal 
minimum number of components present in the full state if we are given the freedom to 
choose a basis from each subspace. If we are given such a choice in the full Hilbert space 
we would have just one component with the basis having one of its directions aligned along 
the eigenstate, and the von Neumann entropy of the pure state is zero. This may also be 
to compared to a M dimensional random matrix eigenvector, belonging to the Gaussian 
Orthogonal ensemble (GOE) [JTJJ whose Shannon entropy is approximately log(0.5M). 

We turn now to a somewhat more detailed study of the entanglement inherent in indi- 
vidual eigenstates. In Fig. @ is plotted the individual entropies corresponding to all the 
states of a fairly chaotic system. While most of the states have already achieved the entropy 
corresponding to the saturation value of log(.59iV) (small dots), there are many states that 
are prominently low in entanglement (those with an entropy less than 3 is marked with a 
circled dot). The Shannon entropy, not displayed here, of the (small dot) states is to a large 
accuracy log(0.5iV 2 ), (as the dimensionality is M = N 2 ), while there are also minima that 
largely match with those in Fig. (|2]) and therefore these are expected to be localized states. 
In Fig. (£J) is shown two wavef unctions, one a typical chaotic state and the other a localized 
state corresponding to the first prominent minima of the entanglement entropy shown in 
Fig. @) (i = 212). 

Many of the low entropic states are similar to the localized state in Fig. (|3|), and suggests 
a "scarring" Jl5| mechanism. In fact this results from a scarring due to a fixed point of 
mixed stability, and may be called semi-scarred if one is to strictly define scarring as due to 
unstable, hyperbolic orbits |L6|]. The initial condition (qi = 0.5, pi = 0,q2 = 0,p2 = 0) is 
a fixed point for all values of the parameters. When 6 = 2 the eigenvalues of the Jacobian 
at this point has a real pair, corresponding to hyperbolic motions and a complex conjugate 
pair corresponding to stable motions. This fixed point is poised to completely lose stability 
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just after 2 (at pa 2.01). While providing a new example of scarring in higher dimensional 
systems due to orbits of mixed stability type, this shows that entanglement is sensitive to 
eigenfunction localization or scarring. We might expect this to be true for a large class of 
strongly scarred states. 

The structure of the reduced density matrices corresponding to these two states are 
displayed in Fig. @ and reveal the connection to decoherence phenomena. The reduced 
density matrices are typically diagonally dominant in the presence of large scale chaos. 
It is interesting that the localized states are more "cleanly" diagonal than the typically 
delocalized state. There is a rather rapid transition from a non-diagonal to a predominantly 
diagonal density matrix as the system undergoes a transition to chaos; even mixed phase 
spaces seem to lead to sufficiently diagonal density matrices. This is shown in Fig. (||) where 
the average (of the absolute value squared) of the density matrices corresponding to the 
entire spectrum is shown. The reduced density matrices of non-stationary states also tend 
to a diagonal structure very rapidly when the system is chaotic, as will be demonstrated 
later. 

Due to the reduced density matrices being diagonally dominant, it is not obvious if 
they themselves have random matrix properties. However the state in the full Hilbert space 
possesses properties of eigenvectors of random matrices of the GOE (or COE) type. We may 
then use this to understand the typical diagonal nature of the reduced density matrices. The 
diagonal term is (dropping the state index i in favor of the indices indicating the position 
in the reduced density matrix): 

Plnn = £ \(nk\4>)\ 2 - (8) 
k 

Typical (ensemble averaged) value of | (n k\ip) | 2 is 1/N 2 , which one can see easily from merely 
the normalization condition. Therefore the typical diagonal element is of order 1/N. The 
"strength" of the diagonal elements as measured in the spectral average, and shown in 
Fig. © will be 



<|Plnn| 2 >~^ (9) 



The off-diagonal element is given by 



Pmn' = "£(nk\i))(i/)\n'k). (10) 

k 



The ensemble average of this quantity vanishes. The strength of these elements is 

1 1 

(N 2 ) 2 ~ Jp 



< |/W| 2 >~<E \(nk\i>)^\n>k)\ 2 >~ Nj^- 2 = ± (11) 

k 



where we have used the GOE result |L4| that 



«\ mW )\>» =w * TT) „±, (MO, (12) 

and M = N 2 is the dimensionality of the matrix. Thus the average off-diagonal element will 
be smaller than the diagonal by a factor of \J~N . This is borne out to a large extent by the 
numerical results. 
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Is there any advantage in using the reduced density matrix, rather than say the Shannon 
entropy as a measure of eigenfunction properties, in particular of localization? Clearly the 
reduced density matrix contains much more by way of information, and the von Neumann 
entropy derived from it is just one piece of information that gives a global idea of localization. 
If we have the complementary reduced density matrix obtained by tracing out the first degree 
of freedom we will have complete information about the wavefunction, via the Schmidt 
decomposition. We now study the spectral properties of the reduced density matrices and 
observe how information about localization is stored in a remarkably graded manner. 

The eigenvalues (Xj',j = 1, . . . , N) of the reduced density matrix, assumed to be arranged 
in increasing order, also naturally contain information about localization. As we noted earlier 
these eigenvalues fall off exponentially as seen in Fig. flU). The localized state seems to have 
at least two exponential scales while the generic delocalized state has only one. For the 
latter class of states we find numerically that 



where 7 is an N independent constant. While this is valid for the most significant few of 
the eigenvalues, there is a clear deviation of this for smaller eigenvalues and the exponential 
law seems modulated by a polynomial one that we will not investigate in more detail. 

The role of the "temperature" is played by the Hilbert space dimensionality N. This may 
be compared to the "temperature of the eigenvalue gas", which measures the equilibrium 
(fully chaotic) distribution of the level velocities, which is also proportional to N, but the 
proportionality constant (the "Boltzmann constant") may be a measure of system specific 
classical correlations ||17|| . The localized state seems to be cleanly split into two parts, one 
with a localization dominant part and the other which behaves like a generic extended state. 
The degree to which there are two distinguishable scales depends on the particular localized 
state. Thus the information about the state's localization is present in the first few Schmidt 
states of the reduced density matrix. 

We now study the corresponding eigenvectors of the reduced density matrices. These 
are the Schmidt states \4>j) in the expansion of Eq. (Q) for eigenstates. Since these are in 
N dimensional subspaces spanned by either degree of freedom, they correspond classically 
to "projections in the (qi,Pi) space". Here % = 1,2 while tracing out the second or the first 
degree of freedom respectively. We see them in two ways: one is the usual position basis, the 
other is a phase-space representation, such as a Husimi distribution. The latter will reveal 
phase-space scarring effects, provided of course that we know enough about the classical 
dynamics. 

The position basis representation of some of the principal eigenvectors or the Schmidt 
vectors of the reduced density matrices of two states is shown in Fig. (^).One of these states 
is a typical non-localized state while the other is the localized state discussed above. While 
the Schmidt vectors of the non-localized state are essentially unremarkable and appear to be 
random, those that belong to the localized state are themselves localized and appear arranged 
in the manner typical of eigenstates with increasing node numbers. Beyond the six states 
shown here, Schmidt vectors of the localized state also look random and indistinguishable 
from the delocalized state, as is already evidenced in the sixth state. This agrees with the 
fact that the eigenvalues of the Schmidt vectors too seem to share common trends beyond 




(13) 
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this point. Thus information about state localization is stored in the Schmidt vectors in a 
graded manner. 

The Husimi of the Schmidt vectors will reveal more about the classical structures that 
influence the localization. We have asserted previously that the localized state we have been 
analyzing is influence by the fixed point (0.5, 0, 0, 0) and while this is plausible from the state 
vector in the position basis already displayed, it is confirmed by the Husimi of the Schmidt 
vectors. Thus we plot 



W{q,p) = \{q,p\<f>) 



S\2 



in Fig. (H) where \q,p) is a coherent state on the two-torus as developed in |L8]. We have 
also plotted the complementary Schmidt vectors from tracing out the first degree of freedom 
and will therefore give the (<?2,P2) coordinates of any classical structure. The fixed point 
(0.5, 0, 0, 0) is clearly seen projected onto the two subspaces in the Husimis, thus confirming 
our earlier statement that the state is "scarred" by this orbit. In general we may expect 
that periodic orbits scarring the states will be seen in their projections in the Husimis of 
their principal Schmidt vectors. 

The full Husimi distribution is of course four-dimensional and taxes our visualization 
abilities. The above may be compared to a similar approach that has already been in use 
when four dimensional Husimis of two-dimensional eigenfunctions of chaotic oscillators with 
two-degrees of freedom were analyzed via their "quantum surface of sections" [J19||. The 



difference is that the Schmidt states provide a much more complete and systematic way of 
analyzing higher dimensional wavefunctions than the somewhat ad-hoc constructions thus 
far in use. It has been shown for two-dimensional maps and two-degree of freedom flows 
(which are equivalent) that the zeros of the Husimi distribution, dubbed as "stellar represen- 
tations", provide a unique description of the eigenfunctions PD[ . The stellar representation 
of Schmidt vectors may provide a way of avoiding complex functions in two variables; as 
each vector and its correlated partner will have N zeros each and there are N such pairs, we 
would have a total of 2 iV 2 zeros. Of course these in themselves are not sufficient to specify 
the state (as we need the eigenvalues of the reduced matrix as well), but we have seen that 
important information about the states is already present in the Schmidt vectors and must 
be reflected in their zeros. Calulations not presented here indeed confirm this. 



B. Time Dependent Properties 

We dwell briefly on time dependent properties. While for wavefunctions on the full 
Hilbert space, time dependent properties may largely be derived from the stationary one, 
the situation is not entirely obvious when we restrict ourselves to reduced density matrices, 
or "shadows" in restricted subspaces. If we take a arbitrary (nonstationary) state |0o) and 
evolve it according to \<f>(T)) = U T \(po). The reduced density matrix at any time T cannot 
be derived based solely on the reduced density matrices of individual states as: 

Pl (T) = tr 2 (|0(T))(0(T)|) = £l<^o>| 2 Pi* + 

k 

^exp(i(Vfc-V'O^)(V'it|0o}(0o|V'/)tr2(|^}(^|) ) (14) 
k+i 
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where ipk are eigenangles of the corresponding eigenstates. If we assume a non-degenerate 
spectrum, the reduced density matrix averaged over all time is the first sum of the above 
expression and is simply a weighted sum over the reduced density matrices of individual 
eigenstates. 

Thus we expect, based on our previous discussion of diagonally dominant density matrices 
of eigenstates, that any arbitrary state's reduced density matrix will rapidly evolve to a 
predominantly diagonal one. This is of course reminiscent of "decoherence" phenomena 
PJ21[] and indeed as far as each degree of freedom individually is concerned, it is an open 
system from which phase information can flow out or decohere. However, the state will evolve 
in such a way that even the diagonal part of the density matrix is significantly altered, i.e., 
it is not comparable to a situation wherein a pure state is "reduced" to a classical ensemble 
by a measurement like process. Fig. ([5]) shows the evolution of an initial state which is 
the sum of two well separated Gaussian in position representation. The initial reduced 
density matrix has large off-diagonal parts. The evolution under a weak coupling produces 
interesting structures. After even one time step, for low couplings we see the density matrix 
having the off-diagonal parts significantly reduced. In fact this picture "looks" like the one 
resulting from decoherence [21]|. However we are doing a numerically exact computation and 
not using any approximate master equation. For later times the case of small coupling or 
low chaos leads to fairly non-diagonal density matrices, while in comparison the right panel 
shows the evolution of the same initial state for the first four time units and one sees that 
effective diagonality is rapidly achieved. 

The time evolving density matrices represent the "Schmidt paths" If we start with 



an initially disentangled state, the dynamics when completely chaotic, will quickly rotate the 
state into those in which there is maximal entanglement in some sense. Is this entanglement 
different from that observed for stationary states? In fact from numerical calculations not 
shown here it seems that these are identical, modulo fluctuations. This is unlike the case 
of the Shannon entropy which for a time evolving pure state is different from a stationary 
state, basically due to the fact that the time evolving state is in general complex while (for 
time-reversal symmetric systems as we are currently discussing) the stationary state is real. 
The typical behaviour of the entropy is as expected a rise from zero to the constant value 
of log(0.59iV) with a rapidity that is a function of the coupling strengths and hence of the 
chaos in the system. 

We briefly also comment on the case when time reversal invariance is broken. That is 
achieved in the model studied in this paper by choosing a non-zero /3, which is equivalent 
to an introduction of a magnetic flux line. While of course much of what has been said 
already carries over to this case, there may be quantitative differences. For one, the average 
entanglement of eigenstates is slightly higher at about log(0.61iV). Thus time-reversal break- 
ing interactions may on the average produce more entanglement, and the reduced density 
matrices are sensitive to time-reversal breaking. 



IV. DISCUSSION 

We have presented a variety of essentially numerical results concerning reduced density 
matrices of chaotic systems. We have seen that entanglement or the von Neumann entropy 
can be a good measure of localization and that the Schmidt states provide graded information 
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about the nature of localization. We have seen how chaos aids entanglement and it is a small 
step to extend the universality observed from quantum chaos to reduced density matrices. 
However we have not sufficiently explored the density matrices to be able to comment on 
their randomness or type of universality we may expect. 

There are various ways in which the analysis can be extended. One important direc- 
tion would be to have more than two coupled maps, or in general many particle interacting 
systems. Such higher dimensional systems force us to face several problems. One is that 
Schmidt decomposition is no longer possible, and entanglement becomes much harder to 
quantify. Ways to measure entanglement in multipartite systems have been proposed re- 
cently that may prove to be useful. The other is the increasing complexity of the classical 
system (if there is one), and the exponentially growing numerical task. In perhaps a related 
vein entanglement inherent in the ground state of a anti-ferromagnetic model has been re- 
cently studied Loss of quantum coherence was studied via simple models in [12| by 



coupling two state systems to larger Hamiltonians with matrix representations having ran- 
domly chosen elements. A natural situation in this context and in the spirit of this paper 
is to look at bipartite quantum chaotic systems as we have but whose dimensionality is 
unequal. 

An important step toward a fundamental understanding of why random matrix modeling 
must be successful at all was achieved by the use of semiclassical methods and periodic orbits 



sums . A further direction would be to derive and use semiclassical orbit sums for "partial 
traces" . This obviously has close links with the Feynman- Vernon path integral treatment of 
quantum dissipation. We note that 

$>xp(-#T)tr 2 (W T ) = - Vfc - 2vrm)p lfc , (15) 

T km 

where ipk is the eigenangle corresponding to the state whose reduced density matrix is p±k- 
The partial trace of the propagator is therefore naturally a quantity of interest. The classical 
orbits that will be the stationary paths will be "partially-periodic" orbits which for a given 
time T connect two configuration space points in subsystem 1 while appearing periodic in 
subsystem 2. Such semiclassical analysis of partial traces may provide deeper understanding 
of various aspects of quantum open systems in general. 
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FIGURES 



FIG. 1. Average entropy as a function of the coupling b. From top to bottom the cases corre- 
spond to N = 25, 20 and 15 respectively. The inset is a log-log plot of the same. 

FIG. 2. Entropy of all the 1600 states for the case K x = 0.1, K 2 = 0.15, b = 2 and N = 40.The 
straight line is at the value log(0.59iV), while the states below an entropy of 3 are indicated by 
circles. The localized states (i = 212) is the second encircled point. 

FIG. 3. Eigenfunctions of a typical state (top, i = 6) and a localized state (bottom, i = 212). 
Parameter values are same as that of the previous figure 

FIG. 4. The reduced density matrices, obtained by tracing out the second degree of freedom, for 
the delocalized state (top, i = 6) and the localized state (bottom, i = 212) of the previous figure. 

FIG. 5. The spectral average of the square of the elements of the reduced density matrices for 
the cases of b = 0.05 (above) and 6 = 2 (below). N = 20, and K\ = 0.1 and K 2 = 0.15 in this 
figure 

FIG. 6. Principal eigenvalues of the two reduced density matrices in Fig. (Q) are shown with the 
open circles and the closed circles corresponding to the delocalized and localized states respectively. 
Note that the scale is log-linear 

FIG. 7. The six principal Schmidt vectors for a typical chaotic state (left) and for the localized 
state (i = 212) (right). 

FIG. 8. The Husimi representation W(q,p) of the principal Schmidt vectors of the localized 
state (i = 212). Top corresponds to tracing out the second degree of freedom and the bottom the 
first. 

FIG. 9. Time evolution of an initial reduced density matrix shown in the left topmost figure. 
The left panels correspond to the case b = 0.1 for evolution over three time steps while the right 
corresponds to b = 2 for evolution over four time steps. 
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